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ABSTRACT 

We have analyzed data from a multi-site campaign to observe oscillations in the F5 star Procyon. The data consist 
of high-precision velocities that we obtained over more than three weeks with eleven telescopes. A new method for 
adjusting the data weights allows us to suppress the sidelobes in the power spectrum. Stacking the power spectrum 
in a so-called echelle diagram reveals two clear ridges that we identify with even and odd values of the angular degree 
(l — and 2, and I = 1 and 3, respectively). We interpret a strong, narrow peak at 446 /iHz that lies close to the 
I — 1 ridge as a mode with mixed character. We show that the frequencies of the ridge centroids and their separations 
are useful diagnostics for asteroseismology. In particular, variations in the large separation appear to indicate a glitch 
in the sound-speed profile at an acoustic depth of ^lOOOs. We list frequencies for 55 modes extracted from the 
data spanning 20 radial orders, a range comparable to the best solar data, which will provide valuable constraints 
for theoretical models. A preliminary comparison with published models shows that the offset between observed and 
calculated frequencies for the radial modes is very different for Procyon than for the Sun and other cool stars. We find 
the mean lifetime of the modes in Procyon to be 1.29^Q'4g days, which is significantly shorter than the 2-4 days seen 
in the Sun. 

Subject headings: stars: individual (Procyon A) — stars: oscillations 
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1. INTRODUCTION 

Tlie success of helioseismology and tlie promise of 
asteroseismology have motivated numerous efforts to 
measure oscillations in solar-type stars. These began 
with ground-based observations (for recen t reviews see 
iBedding fc Kieldsen|[2007l:IAerts et al.ll2008[ ) and now ex- 
tend to space-based photometry, p articularly with the 
CoRoT and K epler Missions (e.g., iMichel et all 120081 : 
ICTilliland et all 20101. 

We have carried out a multi-site spectroscopic cam- 
paign to measure oscillations in the F5 star Procyon A 
(HR 2943; HD 61421; HIP 37279). We obtained high- 
precision velocity observations over more than three 
weeks with eleven telescopes, with almost continu- 
ous coverage for th e central ten days. In Paper I 
(jArentoft et al.l[2008l ) we described the details of the ob- 
servations and data reduction, measured the mean oscil- 
lation amplitudes, gave a crude estimate for the mode 
lifetime and discussed slow variations in the velocity 
curve that we attributed to rotational modulation of ac- 
tive regions. In this paper we describe the procedure 
used to extract the mode parameters, provide a list of 
oscillation frequencies, and give an improved estimate of 
the mode lifetimes. 

2. PROPERTIES OF SOLAR-LIKE OSCILLATIONS 

We begin with a brief summary of the relevant prop- 
erti es of so lar-like oscillations (for reviews see, for exam- 
ple. IBrown fc Gilliland 199 J: IBedding fc KieldseiJ [20031 : 
lChristenseii-Dalsgaarcil i2004) . 

To a good approximation, in main-sequence stars 
the cyclic frequencies of low-degree p-mode oscillations 
are regularly s paced, follow ing the asymptotic relation 
(|Tassoullll980l: IGoughlll986l) : 



(1) 



Here n (the radial order) and I (the angular degree) are 
integers, Aj/ (the large separation) depends on the sound 
travel time across the whole star, Dq is sensitive to the 
sound speed near the core and e is sensitive to the reflec- 
tion properties of the surface layers. It is conventional 
to define three so-called small frequency separations that 
are sensitive to the sound speed in the core: ^1^02 is the 
spacing between adjacent modes with I = and / — 2 
(for which n will differ by 1); Sui^ is the spacing between 
adjacent modes with / = 1 and / = 3 (ditto); and 5vqi 
is the amount by which / = 1 modes are offset from the 
midpoint of the I — modes on either sideE3 To be ex- 
plicit, for a given radial order, n, these separations are 
defined as follows: 



(5i^02 = Vnfl — Vn-1,2 

= \{Vn.O + J^«+l.o) - t^n,l 
6^13 = Vn,l - 



(2) 
(3) 
(4) 



If the asymptotic relation (equation[T]) were to hold ex- 
actly, it would follow that all of these separations would 
be independent of n and that 5vo2 = 6I)o, ^^^13 = lOI^o 
and 5vqi = 2Dq. In practice, equation ([Ij is only an 

One can also define an equivalent quantity, <5j^io, as the offset 
oi I = modes from the midpoint between the surrounding I = 1 
modes, so that = Vn,o — \{yn-i,i + 



approximation. In the Sun and other stars, theoretical 
models and observations show that Dq and e vary 
somewhat with frequency, and also with /. Consequently, 
the small separations also vary with frequency. 

The mode amplitudes are determined by the excita- 
tion and damping, which are stochastic processes involv- 
ing near-surface convection. We typically observe modes 
over a range of frequencies, which in Procyon is especially 
broad (about 400-1400 ^Hz; see Paper I). The observed 
amplitudes also depend on I via various projection factors 
(see Table 1 of Kjcldscn ct al. 2008a). Note in particular 
that velocity measurements are much more sensitive to 
modes with I = 3 than are intensity measurements. The 
mean mode amplitudes are modified for a given observ- 
ing run by the stochastic nature of the excitation, result- 
ing in considerable scatter of the peak heights about the 
envelope. 

Oscillations in the Sun are long-lived compared to their 
periods, which allows their frequencies to be measured 
very precisely. However, the lifetime is not infinite and 
the damping results in each mode in the power spectrum 
being split into multiple peaks under a Lorentzian profile. 
The FWHM of this Lorentzian, which is referred to as 
the linewidth F, is inversely proportional to the mode 
lifetime: F = 1/{ttt). We follow the usual definition that 
r is the time for the mode amplitude to decay by a factor 
of e. The solar value of t for the strongest modes ranges 
from 2 to 4 days, as a de creasing function of frequency 
(e.g., IChaolin et aD[l99l . 

Procyon is an evolved star, with theoretical mod- 
els showing that it is close to, or just past, the end 
of the main sequen c e (e.g ., iGuenth er &: Demargud 



19931 iBaFban et al.l [19991 r Cliabove r et all 119991: 



Di Mauro & Christensen-Dalsgaard 2001; Kervella et al 
2004; Eggc nb crgcr^ c t al. 2005: Provos t et al. .200^ 
[Bonanno et afl l2007t iGuenther'et all MM . As such, 
its oscillation spectrum may show deviations from the 
regular comb-like structure described by equation ([T]), 
especially at low frequencies. This is because some 
modes, particularly those with 1 = 1, are shifted by 
avoided crossings with gravity modes i n the stellar 
core (also called 'ni ode bumping'; see lOsakil 119751 : 
lAizenman et al.lll977[ ). These so-called 'mixed modes' 
have p-mode character near the surface but g-mode 
character in the deep interior. Some of the theoretical 
models of Procyon cited above indeed predict these 
mixed modes, depending on the evolutionary state of the 
star, and we must keep this in mind when attempting to 
identify oscillation modes in the power spectrum. The 
mixed modes are rich in information because they probe 
the stellar core and are very sensitive to age, but they 
complicate the task of mode identification. 

We should also keep in mind that mixed modes are 
expected to have lower amplitudes and longer lifetimes 
(smaller linewidths) than regular p modes because they 
have larger mode inertias (e.g., iChristensen-DalsgaardI 
l2004f ). Hence, for a data series that is many times longer 
than the lifetime of the pure p modes, a mixed mode may 
appear in the power spectrum as a narrow peak that is 
higher than the others, even though its power (amplitude 
squared) is not especially large. 

Another potential complication is that stellar rota- 
tion causes modes with / > 1 to split into multiplets. 
The peaks of these multiplets are characterized by the 
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azimuthal degree m, which takes on values of m = 
0, ±1,...,±/, with a separation that directly measures 
the rotation rate averaged over the region of the star that 
is sampled by the mode. The measurements are particu- 
larly difficult because a long time series is needed to re- 
solve the rotational splittings. We argue in Appendix [X\ 
that the low value of v sin i observed in Procyon implies 
that rotational splitting of frequencies is not measurable 
in our observations. 

3. WEIGHTING THE TIME SERIES 

The time series of velocity observations was obtained 
over 25 days using 11 telescopes at eight observatories 
and contains just over 22 500 data points. As discussed in 
Paper I, the velocity curve shows slow variations that we 
attribute to a combination of instrumental drifts and ro- 
tational modulation of stellar active regions. We have re- 
moved these slow variations by subtracting all the power 
below 280/iIIz, to prevent spectral leakage into higher 
frequencies that would degrade the oscillation spectrum. 
We take this high-pass-filtered time series of velocities, 
together with their associated measurement uncertain- 
ties, as the starting point in our analysis. 

3.1. Noise-optimized weights 

Using weights when analyzing g round-based observa- 
tions of stellar osci llations (e.g., iGilliland et al.l 119931 : 
iFrandsen et al.|[T995l ) allows one to take into account the 
significant variations in data quality during a typical ob- 
serving campaign, especially when two or more telescopes 
are used. The usual practice, which we followed in Pa- 
per I, is to calculate the weights for a time series from the 
measurement uncertainties, at, according to Wi = l/cf ■ 

These "raw" weights can then be adjusted to minimize 
the noise level in the final power spectrum by identifying 
and revising those uncertainties that are too optimistic, 
and at the same time rescaling the uncertainties to be 
in agreement with the actual noise levels in the data. 
This procedure is described in Paper I and references 
therein. The time series of these noise-optimized weights 
is shown in Figure [Ha. These are the same as those shown 
in Figure Id of Paper I, but this time as weights rather 
than uncertainties. 

The power spectrum of Procyon based on these noise- 
optimized weights is shown in Figure[2]a. This is the same 
as shown in Paper I (lower panel of Figure 6), except 
that the power at low frequencies, which arises from the 
slow variations, has been removed. As described in Pa- 
per I, the noise level above 3 mHz in this noise-optimized 
spectrum is 1.9cms~^ in amplitude. This includes some 
degree of spectral leakage from the oscillations and if we 
high-pass filter the spectrum up to 3 mHz to remove the 
oscillation signal, the noise level drops to 1.5cms~^ in 
amplitude. 

The task of extracting oscillation frequencies from the 
power spectrum is complicated by the presence of struc- 
ture in the spectral window, which are caused by gaps or 
otherwise uneven coverage in the time series. The spec- 
tral window using the noise-optimized weights is shown 
in Figure [Da. Prominent sidelobes at ±11.57 //Hz cor- 
respond to aliasing at one cycle per day. Indeed, the 
prospect of reducing these sidelobes is the main reason 
for acquiring multi-site observations. However, even with 
good coverage the velocity precision varies greatly, both 



for a given telescope during the run and from one tele- 
scope to another (see Figure [TJa). As pointed out in Pa- 
per I, using these measurement uncertainties as weights 
has the effect of increasing the sidelobes in the spectral 
window. We now discuss a technique for addressing this 
issue. 

3.2. Sidelobe- optimized weights 

Adjusting the weights allows one to suppress the side- 
lobe structure; the trade-off is an increase in the noise 
level. This technique is routinely used in radio astron- 
omy when synthesising im ages from interferometers (e.g., 
iHo gbom fc Brouwlll974l ). An extreme case is to set all 
weights to be equal, which is the same as not using 
weights at all. This is certainly not optimal because it 
produces a power spectrum with greatly increased noise 
(by a factor of 2.3) but still having significant sidelobes, 
as can be seen in Figure 6a of Paper I. Adjusting the 
weights on a night- by- night basis in order to minimize 
the sidelobes was used in the analysis of du al-site ob- 
servations of a Cen A (jBedding et al.l I2004D. a Cen B 
(|Kieldsen et al.|[2005[) . and d Hvi (jBedding et aL|[200l . 
For our multi-site Procyon data this is impractical be- 
cause of the large number of (partly overlapping) tele- 
scope nights. We have developed a more general algo- 
rithm for adjusting weights to minimize the sidelobes (H. 
Kjeldsen et al., in prep.). The new method, which is su- 
perior because it does not assume the oscillations are 
coherent over the whole observing run, is based on the 
principle that equal weight is given to all segments of the 
time series. The method produces the cleanest possible 
spectral window in terms of suppressing the sidelobes, 
and we have tested it with good r esults using publishe d 
data for a Cen A and B, and 13 Hyi (|Arentoft et al.ll2010D . 

The new method operates with two timescales, Ti and 
T2. All data segments of length Ti (=2hr, in this case) 
are required to have the same total weight throughout the 
time series, with the relaxing condition that variations 
on time scales longer than T2 (=12 hr) are retained. To 
be explicit, the algorithm works as follows. We adjust 
the weights so that all segments of length Ti have the 
same total weight. That is, for each point Wi in the time 
series of weights, define {Si} to be the set of weights 
in a segment of width Ti centered at that time stamp, 
and divide each Wi by the sum of the weights in {Si}. 
However, this adjustment suffers from edge effects, since 
it gives undue weight to points adjacent to a gap. To 
compensate, we also divide by an asymmetry factor 



R = 1 



ight 



^Icft 



ight 



(5) 



Here, Sioft is the sum of the weights in the segment {Si} 
that have time stamps less than ti , and Slight is the sum 
of the weights in the segment {Si} that have time stamps 
greater than ti. Note that R ranges from 1, for points 
that are symmetrically placed in their T bin, up to 2 for 
points at one edge of a gap. 

Once the above procedure is done for Ti, which is 
the shortest timescale on which we wish to adjust the 
weights, we do it again with T2, which is the longest 
timescale for adjusting the weights. Finally, we divide 
the first set of adjusted weights by the second set, and 
this gives the weights that we adopt (Figure [TJ)) . 
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3.3. The sidelobe- optimized power spectrum 

Figure [5]!) shows the power spectrum of Procyon based 
on the sidclobe-optimized weights. The spectral window 
has improved tremendously (Figure ISE*) , while the noise 
level at high frequencies (above 3 mHz) has increased by 
a factor of 2.0. 

The power spectrum now clearly shows a regular series 
of peaks, which are even more obvious after smoothing 
(Figure HJc) . We see that the large separation of the star 
is about 55 /iHz, c onfirming the value indicated by several 
_2revious studies (iMosser et al.lll998HMartic et al.lll999l 
2004; E ggenberger et aLl 12004^ iRegulo fc Roca Cortli 
2005tiLeccia et al.ll2007l:IGuenther et al.ll2008[ ). The verv 



and has been followe d u p by numerous authors 
(iBenomar et al.l l2009airbl : iGruberbauer et"aLl 



strong peak at 446 /xHz appears to be a candidate for a 
mixed mode, especially given its narrowness (see Sec- 
tion H]). 

Plotting the power spectrum in echelle format using a 
large separation of 56 fiRz (Figure |4|) clearly shows two 
ridges, as expected^ The upper parts are vertical but 
the lower parts are tilted, indicating a change in the large 
separation as a function of frequency. This large amount 
of curvature in the echelle diagram goes a long way to- 
wards explaining the lack of agreement between previous 
studies on the mode frequencies of Procyon (see the list 
of references given in the previous paragraph) . 

The advantage of using the sidelobe-optimized weights 
is demonstrated by Figure [5] This is the same as Figure |4] 
but for the noise-optimized weights and the ridges are no 
longer clearly defined. 

4. IDENTIFICATION OF THE RIDGES 

We know from asymptotic theory (see equation[T]) that 
one of the ridges in the echelle diagram (Figure 0]) cor- 
responds to modes with even degree {1 = and 2) and 
the other to modes with odd degree (^ = 1 and 3). How- 
ever, it is not immediately obvious which is which. We 
also need to keep in mind that the asymptotic relation 
in evolved stars does not hold exactly. We designate 
the two possibilities Scenario A, in which the left-hand 
ridge in Figure 0] corresponds to modes with odd degree, 
and Scenario B, in which the same ridge corresponds to 
modes with even degree. Figure IHl shows the Procyon 
power spectrum collapsed along several orders. We see 
now double peaks that suggest the identifications shown, 
which corresponds to Scenario B. 

We can check that the small separation Svqi has the 
expected sign. According to asymptotic theory (equa- 
tion [T|), each / = 1 mode should be at a slightly lower fre- 
quency than the mid-point of the adjacent I = modes. 
This is indeed the case for the identifications given in 
Figure ini but would not be if the even and odd degrees 
were reversed. We should be careful, however, since Siygi 
has b een observed to have the opposite sign in red giant 
stars ([Carrier et al.llMoH Bedding e t al.|[20Tol) . 

The problem of ridge identifi cation in F stars 
was first encountered by lAppourchaux et"aLl ()2008f ) 
when analysing CoRoT observations of HD 49933 

When making an echelle diagram, it is common to plot u ver- 
sus (u mod Ai^), in which case each order slopes upwards slightly. 
However, for gray-scale images it is preferable to keep the orders 
horizontal, as was done in the original presentation of the diagram 
IGrec et al.lll983 ). We have followed that approach in this paper, 
and the value given on the vertical axis indicates the frequency at 
the middle of each order. 



2009; 



Mosser fc Appourc hauxl l2009t iRoxburghl 12009 : 
Kallinger et al. 20l3)- Two other F stars observed 
by CoRoT ha ve presented t he sam e problem, namely 
HP 181906 ([Garcia et al.l [2001 and HD 181420 
([Barban et al.l I2009D. A discussion o f the issue was 
recently given bv IBedding fc KieldsenI ([2010l ). who pro- 
posed a solution to the problem that involves comparing 
two (or more) stars on a single echelle diagram after 
first scaling their frequencies. 

Figure [7| shows the echelle diagram for Procyon 
overlaid with scaled frequencies for two stars ob- 
served by CoRoT, using the method described by 
IBedding fc KieldsenI ([2010 ). The filled symbols are 
scaled oscillation fre quencies for the GO st ar HD 49385 
observed by CoRoT ([Deheuvels et al.ll2010D . The scaling 
involved multiplying all frequencies by a factor of 0.993 
before plotting them, with this factor being chosen to 
align the symbols as closely as possible with the Procyon 
ridges. For this star the CoRoT data gave an unam- 
biguous mode identification, which is indicated by the 
symbol shapes. This confirms that the left-hand ridge of 
Procyon corresponds to modes with even I (Scenario B). 

The open symbols in Figure [7| are oscillation fre- 
quencies for HD 49933 from the revised identifica- 
tion bv IBenomar et ahl ([2009bl Scenario B), after mul- 
tiplying by a scaling factor of 0.6565. The align- 
ment with HD 49385 wa s already demonstrated by 
IBedding k KieldsenI ([Ml . We show HD 49933 here 
for comparison and to draw attention to the different 
amounts of bending at the bottom of the right-hand 
(/ = 1) ridge for the three stars. The CoRoT target 
that is most similar to Procyon is HD170987 but un- 
fortunately the S/N ratio is too low to provid e a clear 
identification of the ridges ([Mathur et al.llMl . 

The above considerations give us confidence that Sce- 
nario B in Procyon is the correct identification, and we 
now proceed on that basis. 

5. FREQUENCIES OF THE RIDGE CENTROIDS 

Our next step in the analysis was to measure the cen- 
troids of the two ridges in the echelle diagram. We 
first removed the strong peak at 446 fjHz (it was re- 
placed by the mean noise level). We believe this to 
be a mixed mode and its extreme power means that it 
would significantly distort the result. We then smoothed 
the power spectrum to a resolution of 10/iHz (FWHM). 
To further improve the visibility of the ridges, we also 
averaged across several orders, which corresponds to 
smoot hing in the vertical direction in the echelle dia- 
gram ([Bedding et al.l 120041 : [Kjeldsen et all 120051 : iKarol 
[20f)l . That is, for a given value of Av we define the 
"order-averaged" power-spectrum to be 



OAPS{iy, Aiy) = c.jPS{v + j Av). (6) 

The coefficients Cj are chosen to give a smoothing with 
a FWHM of kAv: 



1 



C-i = c_ 



1 + (2j7fc)2 



(7) 
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We show in Figure [5] the OAPS based on smoothing 
over 4 orders (fc = 4.0), and so we used (cq, . . . , C4) = 
(1,0.8,0.5,0.31,0.2). The OAPS is plotted for three val- 
ues of the large separations (54, 55 and 56 fiHz) and they 
are superimposed. The three curves are hardly distin- 
guishable and we see that the positions of the maxima 
are not sensitive to the value of Av. 

We next calculated a modified version of the OAPS 
in which the value at each frequency is the maximum 
value of the OAPS over a range of large separations (53- 
57 fiHz). This is basically the same as the comb re sponse 



as used previously by several authors ([K i eldse n et al 



19951 iMosser et a"Llll998t iMartic et aT][T999i: ,Leccia et al 



2007fl . The maxima of this function define the centroids 
of the two ridges, which are shown in Figure |9l 

In Figure [TU] we show the full power spectrum of Pro- 
cyon (using sidelobe-optimized weights) collapsed along 
the ridges. This is similar to Figure |6] except that each 
order was shifted before the summation, so as to align 
the ridge peaks (symbols in Figure |9]) and hence remove 
the curvature. This was done separately for both the 
even- and odd-degree ridges, as shown in the two pan- 
els of Figure [TOl The collapsed spectrum clearly shows 
the power corresponding to I = 0-3, as well as the extra 
power from the mixed modes (for this figure, the peak at 
446 /xHz has not been removed). 

In Section [6] below, we use the ridges to guide our 
identification of the individual modes. First, however, 
we show that some asteroseismological inferences can be 
made solely from the ridges themselves. This is explained 
in more detail in Appendix [Bl 

5.1. Large separation of the ridges 

Figure lllb shows the variation with frequency of the 
large separation for each of the two ridges (diamonds 
and triangles) . The smoothing across orders (equation[6]) 
means that the ridge frequencies are correlated from one 
order to the next and so we used simulations to estimate 
uncertainties for the ridge centroids. 

The oscillatory behavior of Av as a function of fre- 
quency seen in Figure [TTI a is pr esumably a si gnature of 
the helium ionization zone (e.s. IGoughl[l99flh . The os- 
cillation is also seen in Figure [TTb . which shows the sec- 
ond differences for the two ri dges, defined as follows (see 
iGouah 1990: Ballot et al.12004 - .Houdek fc Gough,2007') : 



2^n,even — ^n— l.even '^^n, even ~t~ even (8) 



^2^n,even — ^n— l,even 2z^', _ _ 
A2l^n,odd = ^'n-l,odd " 2i^„^odd + ^'n+l,odd- (9) 

The period of the oscillation is ^500 /iHz, which implies 
a glitch at an acoustic depth that is approximately twice 
the in verse of this value (Gough 1990; Houdek & Gough 
|2007| ). namely ~ 1000 s. To determine this more precisely, 
we calculated the power spectrum of the second differ- 
ences for both the odd and even ridges, and measured 
the highest peak. We found the period of the oscillation 
in the second differences to be 508 ± 18 /iHz. Comparing 
this result with theoretical models will be the subject of 
a future paper. 

The dotted lines in Figure [TTIa show the variation of At^ 
with frequency calculated from t he autocorrelation of the 
time s eries usi ng the method of IMosser fc Appourchau^^ 
(|2009l see also lRoxburgh fc Vorontsovll2006[ ). The mixed 
mode at 446 /xHz was first removed and the smoothing 



filter had FWHM equal to 3 times the mean large sepa- 
ration. We see general agreement with the values calcu- 
lated from the ridge separations. Some of the differences 
presumably arise because the autocorrelation analysis of 
the time series averages the large separation over all de- 
grees. 

5.2. Small separation of the ridges 

Using only the centroids of the ridges, we can measure 
a small separation that is useful for asteroseismology. By 
analogy with Svqi (see equation [3]) , we define it as the 
amount by which the odd ridge is offset from the mid- 
point of the two adjacent even ridges, with a positive 
value corresponding to a leftwards shift (as observed in 
the Sun). That is. 



3dd 



n-f l,even 



^n,odd • 



(10) 



Figur e ITTb shows our measurements of this small sepa- 
ration0 It is related in a simple way to the conven- 
tional small separations 5vqi^ Svo2i and dvi3 (see Ap- 
pendix |B] for details) and so, like them, it gives infor- 
mation about the sound speed in the core. Our mea- 
surements of this small separation can be compared with 
theoretic al models using the equations in Appe ndix IB] 
(e.g., see iChristensen-Dalsgaard fc Houdekl[2009h . 

6. FREQUENCIES OF INDIVIDUAL MODES 

We have extracted oscillation frequencies from the time 
series using the standard procedure of iterative sine- wave 
fitting. Each step of the iteration involves finding the 
strongest peak in the sidelobe-optimized power spectrum 
and subtracting the corresponding sinusoid from the time 
series. Figure [12] shows the result. The two ridges are 
clearly visible but the finite mode lifetime causes many 
modes to be split into two or more peaks. We might 
also be tempted to propose that some of the multiple 
peaks are due to rotational splitting but, as shown in 
Appendix [XI this is unlikely to be the case. 

Deciding on a final list of mode frequencies with cor- 
rect I identifications is somewhat subjective. To guide 
this process, we used the ridge centroids shown in Fig- 
ure [9| as well as the small separations Si'02 and 61^13 from 
the collapsed power spectrum (see Figures [6| and [T0| . 
Each frequency extracted using iterative sine-wave fit- 
ting that lay close to a ridge was assigned an I value and 
multiple peaks from the same mode were averaged. The 
final mode frequencies are listed in Table [1] while peaks 
with S/N > 3.5 that we have not identified are listed in 
Table [21 Figures [T3l and [TTI show these peaks overlaid on 
the sidelobe-optimized power spectrum. Figure [TS] shows 
the three small separations (equations [2H11) as calculated 
from the frequencies listed in Table [T] The uncertain- 
ties in the mode frequencies are shown in parentheses in 
Table [H These depend on the S/N ratio of the peak and 
were calibrated using simulations (e.g., see lBedding et al.l 
I2OOI . 

The entries in Table [5] are mostly false peaks due to 
noise and to residuals from the iterative sine-wave fit- 
ting, but may include some genuine modes. To check 

We could also define a small separation Su^^^ even to be the 
amount by which the centroid of the even ridge is offset rightwards 
from the midpoint of the adjacent odd ridges. This gives similar 
results. 
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whether some of them may be daily ahases of each other 
or of genuine modes, we calculated the differences of all 
combinations of frequencies in Tables [T] and [2] The his- 
togram of these pairwise differences was flat in the vicin- 
ity of 11.6 fiHz and showed no excess, confirming that 
daily aliases do not contribute significantly to the list of 
frequencies in the tables. 

We also checked whether the number peaks in Table [2] 
agrees with expectations. We did this by analysing a 
simulated time series that matched the observations in 
terms of oscillations properties (frequencies, amplitudes 
and mode lifetimes), noise level, window function and 
distribution of weights. We extracted peaks from the 
simulated power spectrum using iterative sine-wave fit- 
ting, as before, and found the number of "extra" peaks 
(not coinciding with the oscillation ridges) to be simi- 
lar to that seen in Figure [121 Finally, we remark that 
the peak at 408 /iHz is a candidate for a mixed mode 
with I — 1, given that it lies in the same order as the 
previously identified mixed mode at 446 fiHz (note that 
we expect one extra I = 1 mode to occur at an avoided 
crossing). 

The modes listed in Table [T] span 20 radial orders and 
more than a factor of 4 in frequency. This range is similar 
to that obtained from long-term studies of the Sun (e.g., 
iBroomhall et al.l I2009D and is unprecedented in astero- 
seismology. It was made possible by the unusually broad 
range of excited modes in Procyon and the high S/N of 
our data. Since the stellar background at low frequencies 
in intensity measurements is expected to be much higher 
than for velocity measurements, it seems unlikely that 
even the best data from the Kepler Mission will return 
such a wide range of frequencies in a single target. 

7. MODE LIFETIMES 

As discussed in Section [21 if the time series is suf- 
ficiently long then damping causes each mode in the 
power spectrum to be split into a series of peaks un- 
der a Lorentzian envelope having FWHM F = 1/{ttt), 
where r is the mode lifetime. Our observations of Pro- 
cyon are not long enough to resolve the modes into clear 
Lorentzians, and instead we see each mode as a small 
number of peaks (sometimes one). Furthermore, the cen- 
troid of these peaks may be offset from th e position of the 
true m ode, as illustrated in Figure 1 of [Anderson et ahl 
(|1990f ). This last feature allows one to use the scatter 
of the extracted frequencies about smooth ridges in the 
echelle diagram, cali brated using simulati o ns, to estimate 
the m ode lifetime ([Kieldsen et al.l [20051 : [Bedding et al.l 
[2007[ ). That method cannot be applied to Procyon be- 
cause the I — and / = 2 ridges are not well-resolved 
and the I — 1 ridge is affected by mixed modes. 

Rather than looking at frequency shifts, we have es- 
timated the mode lifetime from the variations in mode 
amplitudes (again calibrated using simulations). This 
method is less precise but has the advantag e of being in- 
dependent of the niode identifications (e.g.. [Leccia et al.l 
[20071 [Carrier et all [20071: IBedding et al.l 120071) . InTa^ 
per I we calculated the smoothed amplitude curve for 
Procyon in ten 2-day segments and used the fluctua- 
tions about the mean to make a rough estimate of the 
mode lifetime: r = l.StJJdays. We have attempted 
to improve on that estimate by considering the ampli- 
tude fluctuations of individual modes, as has been done 



for th e Sun (e.g..lToutain fc Fr6hlic hlll992[:lBa udin et al.l 
119961: [Chang fc Goughlll998D . but were not able to pro- 
duce well-calibrated results for Procyon. 

Instead, we have measured the "pea kiness" of the 
power spectrum (see [Bedding et aI1l20 07) by calculating 
the ratio between the square of the mean amplitude of 
the 15 highest peaks in the range 500-1300 //Hz (found 
by iterative sine- wave fitting) and the mean power in the 
same frequency range. The value for this ratio from our 
observations of Procyon is 6.9. We made a large number 
of simulations (3600) having a range of mode lifetimes 
and with the observed frequency spectrum, noise level, 
window function and weights. Comparing the simula- 
tions with the observations led to a mode lifetime for 
Procyon of 1.29^^^ days. 

This agrees with the value found in Paper I but is more 
precise, confirming that modes in Procyon are signifi- 
cantly more short-lived than those of the Sun. As dis- 
cussed in Section [21 the d ominant modes in th e Sun have 
hfetimes of 2-4 days fe.g.. [Chaplin et al.[[T997l ). The ten- 
dency for hotter stars to h ave shorter mode lifet imes has 
recently been discussed by [Chaplin et al.l ([2009D . 

8. FITTING TO THE POWER SPECTRUM 

Extracting mode parameters by fitting directly to the 
power spectrum is widely used in helioseismology, where 
the time series extends continuously for months or even 
years, and so the indiv idual modes are well- resolved (e.g., 
[Anderson et al.l [19901) . Mode fitting has not been ap- 
plied to ground-based observations of solar-type oscil- 
lations because these data typically have shorter du- 
rations and significant gaps. Global fitting has been 
carried out on spacecraft data, beginning with the 50- 
d tim e series of a Cen A taken with the WIRE space- 
craft ([Fletcher et al.l [2006D and the 60-d light cu rve of 
HP 49933 from CoRoT (|Appourchaux et al.ll2008l) . Our 
observations of Procyon are much shorter than either of 
these cases but, given the quality of the data and the 
spectral window, we considered it worthwhile to attempt 
a fit. 

Global fits to the Procyon power spectrum were made 
by several of us. Here, we present results from a fit using 
a Bayesian approach (e.g. . IGregorvll2005l ) . which allowed 
us to include in a straightforward way our prior knowl- 
edge of the oscillation properties. The parameters to be 
extracted were the frequencies, heights and linewidths of 
the modes. To obtain the marginal probability distri- 
butions of these parameters and their associated uncer- 
tainties, we employed an APT MCMC (Automated Par- 
allel Tempering Markov Chain Monte Carlo) algorithm. 
It implements the Metropolis-Hastings sampler by per- 
forming a random walk in parameter spac e while d rawing 
samples from the posterior distribution ([Gregorvl [2005). 
Further details of our implementation of the algorithm 
will be given elsewhere (T.L. Campante et al., in prep.). 

The details of the fitting are as follows: 

• The fitting was performed over 17 orders (5-21) 
using the sidelobe-optimized power spectrum. In 
each order we fitted modes with ^ = 0, 1, and 
2, with each individual profile being described 
by a symmetric Lorentzian with FWHM F and 
height H. The mode frequencies were constrained 
to lie close to the ridges and to have only small 
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jumps from one order to the next (a Gaussian prior 
with a — 3/iHz). The S/N ratios of modes with 
I = 3 were too low to permit a fit. In order to take 
their power into account, we included them in the 
model with their frequencies fixed by the asymp- 
totic relation (equation [T]) . 

The data are not good enough to provide a use- 
ful estimate of the linewidth of every mode, or 
even of every order. Therefore, the linewidth was 
parametrized as a linear function of frequency, de- 
fined by two parameters Fgoo and ri2oo, which are 
the values at 600 and 1200 /iHz. These parame- 
ters were determined by the fit, in which both were 
assigned a uniform prior in the range 0-10 /iHz. 

The height of each mode is r elated to the linewidt h 
and amplitude according to (jChaplin et al.l[2005[ ): 



H = 



2^2 
ttT 



(11) 



The amplitudes A of the modes were determined 
as follows. For the radial modes {I = 0) we used 
the smoothed amplitude curve measured from our 
observations, as shown in Figure 10 of Paper I. The 
amplitudes of the non-radial modes {I = 1-3) were 
then calculated from the radi al modes using th e 
ratios given in Table 1 of Kicl dsen et al.l (|2008aD . 
namely So : Si : S2 ■ ^ 1.00 : 1.35 : 1.02 : 0.47. 

• The background was fitted as a flat function. 

• We calculated the rotationally-split profiles of the 
non-radial mode s usin g the description given by 
iGizon fc Solankil (12003). The inclination angle of 
the rotation axis was fixed at 31°, which is the in- 
clination of the binary orbit (jG irard "eFall [2OOOI) 
and, as discussed in Paper I (Section 4.1), is con- 
sistent with the rotational modulation of the ve- 
locity curve. The rotational splitting was fixed at 
0.7 /^Hz, which was chosen to match the observed 
value of vsvui = 3.16kms~^ (lAUende Prieto et al.l 
[2002I) . given the known radius of the star. As dis- 
cussed in Appendix 13 choosing different values for 
the inclination (and hence the splitting) does not 
affect the mode profile, assuming reasonable values 
of the linewidth. 

We carried out the global fit using both scenarios dis- 
cussed in Section m The fit for Scenario B is shown as 
the smooth curve in Figure [13] and the fitted frequen- 
cies are given in Table [3] Note that the mixed mode 
at 446 /xHz was not properly fitted because it lies too far 
from the ridge (see the first bullet point above) . To check 
the agreement with the results discussed in Section [51 we 
examined the differences betweens the frequencies in Ta- 
bles [T] and [3] We found a reduced of 0.74, which 
indicates good agreement. A value less than 1 is not sur- 
prising given that both methods were constrained to find 
modes close to the ridges. 

The fitted linewidths (assumed to be a linear function 
of frequency, as described above) gave mode lifetimes of 
1.5 ±0.4 days at 600 /iHz and 0.6 ±0.3 days at 1200 /iHz. 
These agree with the single value of 1.29lQ'4g days found 
above (Section[7]) , and indicate that the lifetime increases 



towards lower frequencies, as is the case for the Sun and 
for the F-type CoRoT t argets HP 49933 (iBenomar et all 
I2009bf) and HD 181420 (|Barban et al.ll2009f) . 

We also carried out the global fit using Scenario A. We 
found through Bayesian model selection that Scenario A 
was statistically favored over Scenario B by a factor of 
10:1. This fact or classifies as "significant" on the scale of 
IJeffrevd (fl96it see Table 1 of lLiddldl2009l) . On the same 
scale, posterior odds of at least ~13:1 are required for a 
classification of "strong to very strong", and "decisive" 
requires at least ^150:1. In our Bayesian fit to Procyon, 
the odds ratio in favor of Scenario A did not exceed 13;1, 
even when different sets of priors were imposed. 

In light of the strong arguments given in Section |4] in 
favour of Scenario B, we do not consider the result from 
Bayesian model selection to be sufficiently compelling to 
cause us to reverse our identification. Of course, it is 
possible that Scenario A is correct and, for complete- 
ness, we show these fitted frequencies in Tabled The fit 
using Scenario A gave mode lifetimes of 0.9 ± 0.2 days at 
600 /iHz and 1.0 ± 0.3 days at 1200 /iHz. 

9. PRELIMINARY COMPARISON WITH MODELS 

A detailed comparison of the observed frequencies of 
Procyon with theoretical models is beyond the scope of 
this paper, but we will make some preliminary com- 
ments on the systematic offset between the two. It is 
well-established that incorrect modeling of the surface 
layers of the Sun is responsible for discrepancies be- 
tween the obser ved and calculated o s cillation frequencies 
{ Christcnscn-D alsgaard et aLl 119881: Dziembowski et all 
il988i: .Rosenthal et al.lll999l: iLi et al.ll2002|)^ 



To ad dress this problem for other stars. iKieldsen etahl 
(|2008bf ) proposed an empirical correction to be applied 
to model frequencies that takes advantage of the fact 
that the offset between observations and models is inde- 
pendent of I and goes to zero with decreasing frequency. 
They measured the offset for the Sun to be a power law 
with exponent b = 4.9 and applied this correction to 
the radial modes of other stars, finding very good results 
that allowed them to estimate mean stellar densities very 
accurately (better than 0.5 per cent). 

We have applied this method to Procyon, comparing 
our observed frequencies for the radial modes with var- 
ious published mod els to determine the sc aling factor r 
and the offset (see IKieldsen et al.l I2008bl for details of 
the method). The results are shown in Figure [TOl Inter- 
estingly, the offset between the observations and scaled 
models does not go to zero with decreasing frequency. 
T his contrasts with the G and K-type stars investigated 
bv lKieldsen et al.l ()2008b[ ). namely the Sun, a Cen A and 
B, and ;9 Hyi. 

The method of'K ieldsen etall ()2008bf ) assumes the cor- 
rection to be applied to the models to have the same form 
as in the Sun, namely a power law with an exponent of 
b — 4.9. The fit in Figure[Tn]is poor and is not improved 
by modest adjustments to h. Instead, the results seem 
to imply an offset that is constant. Setting 6 = and 
repeating the calculations produces the results shown in 
Figure 1171 where we indeed see a roughly constant off- 
set between the models and the observations of about 
20AiHz. 

As a check, we can consider the density implied for Pro- 
cyon. The stellar radius can be calculated from the inter- 
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ferometric radius and t he parallax. The angula r diame- 
ter of 5.404±0.031nias (|Aufdenberg et al.|[2005l Table 7) 
and the rev ised Hipparcos parallax of 285.93 ± 0.88 mas 
(jvan Leeuw en 2007) give a radius of 2.041 ± 0.015 Rq. 

Procyon is in a binary system (the secondary is 
a white dwarf), allowing the mass to be determined 
from astrometry. iGirard et al.l (I2000D fou n d a v alue of 
1.497 ± 0.037 M^. whil e IGatewood fc HmI (l2006h found 
1.431 ± 0.034 Mq (see iGuenther et al.l l2008r for further 
discussion) . 

The density obtained using the fits shown in Figure [16] 
is in the range 0. 255-0. 258gcm~'^. Combining with the 
radius implies a mass in the range 1.54-1.56 M©. The 
density obtained using the fits shown in Figure [17] is in 
the range 0. 242-0. 244gcm~'^, implying a mass of 1.46- 
1.48 Mq. The latter case seems to be in much bet- 
ter agreement with the astrometrically determined mass, 
lending some support to the idea that the offset is con- 
stant. 

We can also consider the possibility that our mode 
identification is wrong and that Scenario A is the correct 
one (see Sections [4] and [8|) . With this reversed identi- 
fication, the radial modes in Procyon are those in Ta- 
ble [1] listed as having I — 1. Assuming these to be radial 
modes, the offset between them and the model frequen- 
cies is again constant, as we would expect, but this time 
with a mean value close to zero. The implied density for 
Procyon is again consistent with the observed mass and 
radius. 

The preceding discussion makes it clear that the cor- 
rection that needs to be applied to models of Procyon is 
very different from that for the Sun and other cool stars, 
regardless of whether Scenario B or A is correct. In par- 
ticular, the substantial nearly-constant offset implied by 
Figure [TBI would indicate errors in the modeling extend- 
ing well beyond the near-surface layers. We also note 
that in terms of the asymptotic expression (equation [T} 
a constant offset would imply an error in the calculation 
of e. 

10. CONCLUSION 

We have analyzed results from a multi-site campaign 
on Procyon that obtained high-pre cision velocity obser- 
vations over more than three weeks (jArentoft et al.|[2008l 
Paper I). We developed a new method for adjusting the 
weights in the time series that allowed us to minimize the 
sidelobes in the power spectrum that arise from diurnal 
gaps and so to construct an echelle diagram that shows 
two clear ridges of power. To identify the odd and even 
ridges, we summed the power across several orders. We 
found structures characteristic of Z = and 2 in one ridge 
and / = 1 and 3 in the other. This identification was con- 
firmed by comparing ou r Procyon da ta in a scaled echelle 
diagram (Bedding fc Kjerdsenll2010l ) with other stars for 
which the ridge identification is known. We showed that 
the frequencies of the ridge centroids and their large and 
small separations are easily measured and are useful diag- 
nostics for asteroseismology. In particular, an oscillation 
in the large separation appears to indicate a glitch in the 
sound-speed profile at an acoustic depth of '--^lOOOs. 

We identify a strong narrow peak at 446 /xHz, which 
falls slightly away from the I = 1 ridge, as a mixed mode. 
In Table [T] we give frequencies, extracted using iterative 
sine-wave fitting, for 55 modes with angular degrees I of 



0-3. These cover 20 radial orders and a factor of more 
than 4 in frequency, which reflects the broad range of 
excited modes in Procyon and the high S/N of our data, 
especially at low frequencies. Intensity measurements 
will suffer from a much higher stellar background at low 
frequencies, making it unlikely that even the best data 
from the Kepler Mission will yield the wide range of fre- 
quencies found here. This is a strong argument in favor of 
continuing efforts towards ground-based Doppler studies, 
such as the SONG network (Stellar Observations Net- 
work Group; Grundahl ct al. 2008), which is currently 
under construction, and the proposed Antarctic instru- 
ment SIAMOIS (Seismic Interfer ometer to Measure Os- 
cillations in the Interior of Stars: [Mosser et aI]|2008D . 

We estimated the mean lifetime of the modes by com- 
paring the "peakiness" of the power spectrum with sim- 
ulations and found a value of 1.29!to 49 days, significantly 
below that of the Sun. A global fit to the power spectrum 
using Bayesian methods confirmed this result and pro- 
vided evidence that the lifetime increases towards lower 
frequencies. It also casts some doubt on the mode iden- 
tifications. We still favor the identification discussed 
above, but leave open the possibility that this may need 
to be reversed. Finally, comparing the observed frequen- 
cies of radial modes in Procyon with published theoreti- 
cal models showed an offset that appears to be constant 
with frequency, making it very different from that seen 
in the Sun and other cool stars. Detailed comparisons of 
our results with theoretical models will be carried out in 
future papers. 

We would be happy to make the data presented in this 
paper available on request. 
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Table 1 

Oscillation Frequencies in Procyon (in /iHz) 



Order 


I = 





I = 


1 


I = 


2 


/ = 


3 


4 






331.3 


(0.8) 










5 






387.7 


(0.7 ) 










6 


415.5 


(0.8) 


445.8 


(0.3) 


411.7 


(0.7) 


■_ 




7 


466.5 


(1.0) 


498.6 


(0.7) 


464.5 


(0.9) 


488.7 


(0.9) 


8 






551.5 


(0.7) 






544.4 


(0.9) 


9 


576.0 


(0.7) 


608.2 


(0.5) 










10 


630.7 


(0.6) 


660.6 


(0.7) 


627.0 


(1.1) 


653.6 


(0.8) 


11 


685.6 


(0.7) 


712.1 


(0.5) 


681.9 


(0.7) 




12 


739.2 


(0.7) 


766.5 


(0.5) 


7,36.2 


Co 5') 






13 


793.7 


(0.9) 


817.2 


(0.6) 


792.3 


(0.9) 






14 


849.1 


(0.7) 


873.5 


(0.6) 


845.4 


(0.6) 


869.5 


(0.6) 


15 


901.9 


(0.8) 


929.2 


(0.7) 




926.6 


(0.6) 


16 


957.8 


(0.6) 


985.3 


(0.7) 


956.4 


(0.5) 


980.4 


(0.9) 


17 


1015.8 


(0.6) 


1040.0 


(0.7) 




1034.5 


(0.7) 


18 


1073.9 


(0.7) 


1096.5 


(0.7) 


1068.5 


(0.7) 




19 


1126.7 


(0.5) 


1154.6 


(0.9) 


1124.3 


(0.9) 






20 


1182.0 


(0.7) 


1208.5 


(0.6) 


1179.9 


(1.0) 






21 


1238.3 


(0.9) 


1264.6 


(1.0) 


1237.0 


(0.8) 






22 


1295.2 


(1.0) 






1292.8 


(1.0) 






23 


1352.6 


(1.1) 


1375.7 


(1.0) 


1348.2 


(1.0) 







Table 2 

Unidentified Peaks 
with S/N > 3.5 







S/N 


(AtHz) 




407.6 


(0.8) 


3.5 


512.8 


(0.8) 


3.6 


622.8 


(0.6) 


4.3 


679.1 


(0.7) 


4.0 


723.5 


(0.6) 


4.7 


770.5 


(0.7) 


4.1 


878.5 


(0.6) 


4.4 


890.8 


(0.7) 


3.6 


935.6 


(0.7) 


3.9 


1057.2 


(0.7) 


3.7 


1384.3 


(0.7) 


3.6 



Table 3 

Frequencies from global fit using Scenario B (in /xHz, with — /+ 
uncertainties) 



Order 


I 


— 







I 


= 1 


I 




2 




5 


363.6 


(0 


8/0 


9) 


387.5 


(0.6/0.6) 


358.5 


(1 


3/1 


2) 


6 


415.3 


(3 


3/1 


0) 






408.1 


(1 


0/3 


7) 


7 


469.7 


(1 


6/2 


1) 


498.8 


(0.7/0.8) 
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Figure 1. Weights for time series of velocity observations of Procyon, optimized to minimize: (a) the noise level and (6) the height of the 
sidelobes. 
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Figure 2. Power spectrum of oscillations in Procyon: (a) using the noise-optimized weights; (6) using the sidelobe-optimized weights; 
(c) using the sidelobe-optimized weights and smoothing by convolution with a Gaussian with FWHM 2 //Hz. 
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Figure 3. Spectral window for tlie Procyon observations using 
(a) noise-optimized weights and (6) sidelobe-optimized weights. 
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Figure 4. Power spectrum of Procyon in echelle format using 
a large separation of 56/.tHz, based on the sidelobe-optimized 
weights. Two ridges are clearly visible. The upper parts are verti- 
cal but the lower parts are tilted, indicating a change in the large 
separation as a function of frequency. The orders are numbered 
sequentially on the right-hand side. 
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Figure 5. Same as Fig. |4l but for the noise-optimized weights. 
The sidelobes from daily aliasing mean that the ridges can no longer 
be clearly distinguished. 
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Figure 6. The power spectrum of Procyon collapsed along several 
orders. Note that the power spectrum was first smoothed slightly 
by convolving with a Gaussian with FWHM 0.5 fiHz. The dotted 
lines are separated by exactly Au/2, to guide the eye in assessing 
the 0—1 small separation 
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Figure 7. Echelle diagram for Procyon smoothed to 2 /xHz 
(greyscale) overlaid with scaled frequencies for two stars ob- 
served by CoRoT. Th e filled sym bols are oscillation frequencies for 
HD 49385 reported by Deheuvels ct al. (2010), after multiplying by 
0.993. Open symbols are oscillation frequencies fo r H D 49933 from 
the revised identification by Benomar et al. ( 2009H, Scenario B) 
after multiplying by 0.6565. Symbol shapes indicate mode degree: 
I = (circles), 1 = 1 (triangles), and I = 2 (squares). 
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Figure 9. Centroids of the two ridges, as measured from the 
comb response. The grayscale shows the sidelobe-optimizcd power 
spectrum from which the peaks were calculated. 
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Figure 8. Order-averaged power spectrum (OAPS), where 
smoothing was done with a FWHM of 4.0 orders (see text). The 
OAPS is plotted for three values of the large separations (54, 55 
and 56 fiRz) and we see that the positions of the maxima are not 
very sensitive to the value of Au. 
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Figure 10. The power spectrum of Procyon collapsed along the 
ridges, over the full range of oscillations (18 orders). The upper 
panel shows the left-hand ridge, which we identify with modes hav- 
ing even degree, and the lower panel shows the right-hand ridge 
(odd degree). Note that the power spectrum was first smoothed 
slightly by convolving with a Gaussian with FWHM 0.6 /iRz. 
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Figure 11. Symbols show the frequency separations in Procyon 
as a function of frequency, as measured from the ridge centroids: 
(a) large frequency separation, (b) second differences, and (c) small 
frequency separation. Tiie dotted lines in panel a show the varia- 
tion in Au (with itlo" range) calculated from the autocorrelation 
of the time series — see the text. 
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Figure 12. Peaks extracted from sidclobe-optimized power spec- 
trum using iterative sine-wave fitting. Symbol size is proportional 
to amplitude (after the ba<;kground noise has been subtracted). 
The grayscale shows the sidelobe-optimized power spectrum on 
which the fitting was performed, to guide the eye. 
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Figure 13. The power spectrum of Procyon at full resolution, with the orders in each column arranged from top to bottom, for easy 
comparison with the echelle diagrams. Vortical dashed lines show the mode frequencies listed in Table 111 and dotted lines show the peaks 
that have not been identified, as listed in Table [2] The smooth curve shows the global fit to the power spectrum for Scenario B (see 
Section (Sjl. 
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Figure 16. The difference between observed frequencies of radial 
modes in Procyon and those of scaled models. The symbols in- 
dicate different models, as fo llows: squares from [ Chaboy er et al.[ 
(1999, Table 2), crosses from 'Di Mauro fc Christensen-DalsgaardI 
1 25511) . a sterisks fr om Kcrvella et al. (2004j Table 4), and trian- 
gles from lEg genb ergcr et al.i (i2005i . model Mia). In each case, the 
dotted cu rve shows the correction calculated using equation (4) of 
IKieldsen et al., (i2008bi) . 



Figure 14. The power spectrum of Procyon overlaid with mode 
frequencies listed in Table [T] Symbols indicate angular degree 
(squares: 1 = 0; diamonds: 1 = 1; crosses: I = 2; pluses: I = 3). 
Asterisks show the peaks that have not been identified, as listed in 
Table [2] 
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Figure 15. Small frequency separations in Procyon, as measured 
from the mode frequencies listed in Table [Tl 



16 



Bedding et al. 



REFERENCES 



Aerts, C, Christensen-Dalsgaard, J., Cunha, M., & Kurtz, D. W. 

2008, Sol. Phys., 251, 3 

Aizenman, M., Smeyers, P., & Weigert, A. 1977, A&A, 58, 41 
AUende Prieto, C, Asplund, M., Lopez, R. J. G., & Lambert, 

D. L. 2002, ApJ, 567, 544 
Anderson, E. R., Duvall, T. L., & Jefferies, S. M. 1990, ApJ, 364, 

699 

Appourchaux, T., et al. 2008, A&A, 488, 705 
Arentoft, T., Kjeldsen, H., & Bedding, T. R. 2010, in GONG 
2008/SOHO XXI Meeting on Solar-Stellar Dynamos asa 
Revealed by Hello- and Asterseismology, ed. M. Dikpati, 
I. Gonzalez-Hernandez, T. Arentoft, & F. Hill (ASP Conf. 
Ser.), in press (arXiv:0901.3632 ) 
Arentoft, T., Kjeldsen, H., Bedding, T. R., et al. 2008, ApJ, 687, 
1180, (Paper I) 

Aufdenberg, J. P., Ludwig, H.-G., & Kervella, P. 2005, ApJ, 633, 
424 

Ballot, J., Garcia, R. A., & Lambert, P. 2006, MNRAS, 369, 1281 
Ballot, J., Turck-Chieze, S., & Garcia, R. A. 2004, A&A, 423, 
1051 

Barban, C., Michel, E., Martic, M., Schmitt, J., Lebrun, J. C., 

Baglin, A., & Bertaux, J. L. 1999, A&A, 350, 617 
Barban, C., et al. 2009, A&A, 506, 51 

Baudin, F., Gabriel, A., Gibert, D., Palle, P. L., & Regulo, C. 

1996, A&A, 311, 1024 
Bedding, T. R., & Kjeldsen, H. 2003, PASA, 20, 203 
Bedding, T. R., & Kjeldsen, H. 2007, in Unsolved Problems in 

Stellar Physics: A Conference in Honour of Douglas Cough, 

AIP Conf. Proc, vol. 948, ed. R. J. Stancliffe, G. Houdek, 

R. G. Martin, & C. A. Tout, 117 
— . 2010, Commun. Asteroseismology, 161, 3 

Bedding, T. R., Kjeldsen, H., Arentoft, T., et al. 2007, ApJ, 663, 
1315 

Bedding, T. R., Kjeldsen, H., Butler, R. P., McCarthy, C., Marcy, 
G. W., O'Toole, S. J., Tinney, C. G., & Wright, J. T. 2004, 
ApJ, 614, 380 

Bedding, T. R., et al. 2010, ApJ Lett., in press l larXiv: 1001. 02291 1 
Benomar, O., Appourchaux, T., & Baudin, F. 2009a, A&A, 506, 
15 

Benomar, O., et al. 2009b, A&A, 507, L13 

Bonanno, A., Kiiker, M., & Paterno, L. 2007, A&A, 462, 1031 
Broomhall, A.-M., Chaplin, W. J., Davies, G. R., Elsworth, Y., 

Fletcher, S. T., Hale, S. J., Miller, B., & New, R. 2009, 

MNRAS, 396, LlOO 
Brown, T. M., & Gilliland, R. L. 1994, ARA&A, 32, 37 
Carrier, F., et al. 2007, A&A, 470, 1059 
— . 2010, A&A, 55, A73 

Chaboyer, B., Demarque, P., & Guenther, D. B. 1999, ApJ, 525, 
L41 

Chang, H.-Y., & Gough, D. O. 1998, Sol. Phys., 181, 251 
Chaplin, W. J., Elsworth, Y., Isaak, G. R., McLeod, C. P., Miller, 

B. A., & New, R. 1997, MNRAS, 288, 623 
Chaplin, W. J., Houdek, G., Elsworth, Y., Gough, D. O., Isaak, 

G. R., & New, R. 2005, MNRAS, 360, 859 
Chaplin, W. J., Houdek, G., Karoff, C, Elsworth, Y., & New, R. 

2009, A&A, 500, L21 
Christensen-Dalsgaard, J. 2004, Sol. Phys., 220, 137 
Christensen-Dalsgaard, J., Dappen, W., & Lebreton, Y. 1988, 

Nat, 336, 634 

Christensen-Dalsgaard, J., & Houdek, G. 2009, Ap&SS, in press 
(arXiv:0911.4629) 

Deheuvels, S., et ah 2010, A&A, submitted 

Di Mauro, M. P., & Christensen-Dalsgaard, J. 2001, in lAU 
Symposium 203: Recent Insights into the Physics of the Sun 
and Heliosphere: Highlights from SOHO and Other Space 
Missions, ed. P. Brekke, B. Fleck, & J. B. Gurman (ASP Conf. 
Ser.), 94 

Dziembowski, W. A., Paterno, L., & Ventura, R. 1988, A&A, 200, 
213 

Eggenberger, P., Carrier, F., & Bouchy, F. 2005, NewA, 10, 195 
Eggenberger, P., Carrier, F., Bouchy, F., & Blecha, A. 2004, 
A&A, 422, 247 



Fletcher, S. T., Chaplin, W. J., Elsworth, Y., Schou, J., & Buzasi, 

D. 2006, MNRAS, 371, 935 
Frandsen, S., Jones, A., Kjeldsen, H., Viskum, M., Hjorth, J., 

Andersen, N. H., & Thomsen, B. 1995, A&A, 301, 123 
Garcia, R. A., et al. 2009, A&A, 506, 41 
Gatewood, G., & Han, I. 2006, AJ, 131, 1015 
Gilliland, R. L., et al. 1993, AJ, 106, 2441 
Gilliland, R. L., et al. 2010, PASP, 122, 131 
Girard, T. M., Wu, H., Lee, J. T., Dyson, S. E., van Altena, 

W. F., et al. 2000, AJ, 119, 2428 
Gizon, L., & Solanki, S. K. 2003, ApJ, 589, 1009 
Gough, D. O. 1986, in Hydrodynamic and Magnetodynamic 

Problems in the Sun and Stars, ed. Y. Osaki (Tokyo: Uni. of 

Tokyo Press), 117 
Gough, D. O. 1990, in Lecture Notes in Physics, Vol. 367, 

Progress of Seismology of the Sun and Stars, ed. Y. Osaki & 

H. Shibahashi (Berlin: Springer), 283 
Grec, G., Fossat, E., & Pomerantz, M. A. 1983, Sol. Phys., 82, 55 
Gregory, P. C. 2005, Bayesian Logical Data Analysis for the 

Physical Sciences (Cambridge University Press) 
Gruberbauer, M., Kallinger, T., Weiss, W. W., & Guenther, D. B. 

2009, A&A, 506, 1043 
Grundahl, F., Christensen-Dalsgaard, J., Arentoft, T., Frandsen, 

S., Kjeldsen, H., j0rgensen, U. G., & Kjaergaard, P. 2008, 

Commun. Asteroseismology, 157, 273 
Guenther, D. B., & Demarque, P. 1993, ApJ, 405, 298 
Guenther, D. B., et al. 2008, ApJ, 687, 1448 
Hogbom, J. A., & Brouw, W. N. 1974, A&A, 33, 289 
Houdek, G., & Gough, D. O. 2007, MNRAS, 375, 861 
Jeffreys, H. 1961, Theory of Probability, 3rd edn. (New York: 

Oxford University Press) 
Kallinger, T., Gruberbauer, M., Guenther, D. B., Fossati, L., & 

Weiss, W. W. 2010, A&A, 510, A106 
Karoff, C. 2007, MNRAS, 381, 1001 

Kervella, P., Thevenin, F., Morel, P., Berthomieu, G., Borde, P., 

& Provost, J. 2004, A&A, 413, 251 
Kjeldsen, H., Bedding, T. R., Arentoft, T., et al. 2008a, ApJ, 682, 

1370 

Kjeldsen, H., Bedding, T. R., & Christensen-Dalsgaard, J. 2008b, 
ApJ, 683, L175 

Kjeldsen, H., Bedding, T. R., Viskum, M., & Frandsen, S. 1995, 

AJ, 109, 1313 
Kjeldsen, H., et al. 2005, ApJ, 635, 1281 

Leccia, S., Kjeldsen, H., Bonanno, A., Claudi, R. U., Ventura, R., 

& Paterno, L. 2007, A&A, 464, 1059 
Li, L. H., Robinson, F. J., Demarque, P., Sofia, S., & Guenther, 

D. B. 2002, ApJ, 567, 1192 
Liddle, A. R. 2009, Annual Review of Nuclear and Particle 

Science, 59, 95 

Martic, M., Lebrun, J.-C, Appourchaux, T., & Korzennik, S. G. 

2004, A&A, 418, 295 
Martic, M., et al. 1999, A&A, 351, 993 
Mathur, S., et al. 2010, A&A, submitted 
Michel, E., et al. 2008, Science, 322, 558 
Mosser, B., & Appourchaux, T. 2009, A&A, 508, 877 
Mosser, B., Appourchaux, T., Catala, C, Buey, J., & SIAMOIS 

team. 2008, Journal of Physics Conference Series, 118, 012042 
Mosser, B., Maillard, J. P., Mekarnia, D., & Gay, J. 1998, A&A, 

340, 457 
Osaki, J. 1975, PASJ, 27, 237 

Provost, J., Berthomieu, G., Martic, M., & Morel, P. 2006, A&A, 
460, 759 

Regulo, C, & Roca Cortes, T. 2005, A&A, 444, L5 

Rosenthal, C. S., Christensen-Dalsgaard, J., Nordlund, A., Stein, 

R. F., & Trampedach, R. 1999, A&A, 351, 689 
Roxburgh, I. W. 2009, A&A, 506, 435 

Roxburgh, I. W., & Vorontsov, S. V. 2006, MNRAS, 369, 1491 
Tassoul, M. 1980, ApJS, 43, 469 
Toutain, T., & Frohlich, C. 1992, A&A, 257, 287 
van Leeuwen, F. 2007, Hipparcos, the New Reduction of the Raw 
Data (Springer: Dordrecht) 



Oscillations in Procyon. II. Frequencies 



17 



APPENDIX 
A. ROTATIONAL SPLITTING 

We expect non-radial modes to be split due to the rotation of the star. The rotation period of Procyon is not known, 
although slow variations in our velocity observations (Paper I) indicated a value of either 1 0.3 days or twice that value. 
The projected rotational velocity has been measured spectroscopically. lAllende Prieto et al. (2002) determined a value 
o f vsini = 3.16 ± . 50 km s~^, although they note that the actual value may be lower by about 0.5 km s^^. 

iGizon fc Solankil ()2003[ ) have studie d the effect of rota tion on the profiles of solar-like oscillations as a function of 
inclination and mode lifetime (see also lBallot et al.ll2006l ). We have repeated their calculations for our observations of 
Procyon (with sidelobe-optimized weights) . The results are shown in Figure [THl which shows the effects of rotational 
splitting, inclination angle and mode lifetime on the theoretical profile of the modesH Note that the calculations do not 
include the stochastic nature of the excitation and so the function shown here should properly be calle d the expectation 
value of the power spectrum, also known as the limit spectrum. Figure ITSl is similar to Figure 2 of IGizon fc Solankil 
(|2003l ) except that instead of fixing the rotation period, we have fixed v sin i to be the measured value. For I = the 
profile does not depend on the inclination angle, while for / = 1, 2 and 3 the solid and dashed lines show calculations 
for i = 30° (Prot — 16.4 days) and i = 80° (Prot = 32.3 days), respectively. In each panel, results are shown for three 
values of the mode lifetime: 1.5 days (top), 3 days (middle) and infinite (bottom). For each mode lifetime, the curves 
for different i and / are all normalized to have the same area. 

We see from Figure [18] that for a fixed usini, the width of the profile stays roughly constant as a function of 
inclination. If the rotation axis of the star happens to be in the plane of the sky (i = 90° ) then the rotation period 
is too low to produce a measurable splitting. At the other extreme, if the inclination is small (so that the rotation is 
close to pole-on), then the rotational splitting will be large but most of the power will be in the central peak (m = 0). 
Either way, once the profile has been broadened by the mode lifetime, the splitting will be unobservable. 

We conclude that for realistic values of the mode lifetime, our observations are not long enough to detect rotational 
splitting in Procyon. The line profiles are broadened by rotation, but it is not possible to disentangle the rotation rate 
from the inclination angle. Rotational splitting is not measurable in Procyon, except perhaps with an extremely long 
data set. The detection of rotational splitting requires choosing a star with a larger wsini or a longer mode lifetime, 
or both. 



B. RELATING RIDGE CENTROIDS TO MODE FREQUENCIES 

As discussed in Section [SJ the frequencies of the ridge centroids are useful for asteroseismology in cases where it is 
difficult to resolve the ridges into their component modes. In this appendix, we relate the frequencies of the ridge 
centroids to those of the underlying modes, which allows us to express the small separation of the ridges (equation ITUl) 
in terms of the conventional small separations (^I'oi, (51^02, and (Ji'is). These relationships will allow the observations 
to be compared with theoretical models. 

The ridge centroids depend on the relative contributions of modes with / = 0, 1, 2, and 3. The power in the even 
ridge is approximately equally divided between I = and 1 — 1, while the odd ridge is dominated by / = 1 but with 
some c ontribution from I = 3. The exact ratios depend on the observing method, as discussed by iKieldsen et al.l 
(l2008aD. For velocity m easurements, such as those presented in this paper for Procyon, the amplitude ratios given by 
IKieldsen et"all (|2008aL their Table 1) yield the following expressions for the centroids in power: 

<evcn =0-49z.„,0 + 0.5lM„_i,2 (Bl) 

<odd = 0.89i/„,i+0.11i/„_i,3, (B2) 

where the superscript indicates these apply to velocity measurements. 

For photometric measurements, such as those currently being obtain ed with the CoRoT aii d Kepler Missions, the 
relative contributions from the various I values are different. Table 1 of IKieldsen eFall H2008aD gives response factors 
for intensity measurements in the three VIRGO passbands, namely 402, 500 and 862 nm. For CoRoT a nd Kepler, it 
is appropriate to use a central wavelength of 650 nm. Using the same method as IKieldsen et al] ()2008al) . we find the 
ratios (in amplitude) for this case to be So '■ Si : S2 '■ S-^ — 1.00 : 1.23 : 0.71 : 0.14. The ridge centroids measured from 
such data would then be 

^''evc„ = 0.66i/„,o + 0.34i/„_i,2 (B3) 
J^n'odd = 0-99z/„,i + 0.0l!.„_i,3, (B4) 

We can express the new small separation of the ridge centroids (equation I10|) in terms of the conventional ones. For 
velocity we have 

5<f.n,oM = '^'^01 - 0.51Ji/02 + 0.11,5i.i3. (B5) 

and for photometry we have 

64len odd = 5^01 - 0.34^1/02 + 0.01,5l^i3. (B6) 



36 



Note that we have made the quite reasonable assumption that the internal rotation has a similar period to the surface rotation. 
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Figure 18. Theoretical line profiles showing rotational splitting for different mode degrees, similar to Fig. 2 of IGizon fc Solankil 1 12003 ) 
but here using a fixed value of vs'mi, namely 3.16kms~^, as measured for Procyon (AUcnde Pricto et al. 2002)). For I = the profile does 
not depend on the inclination angle, while for I = I, 2 and 3 the solid and dashed lines show calculations for i = 30° (Prot = 16.4 days) 
and i = 80° (Prot = 32.3 days), respectively. In each panel, results are shown for three values of the mode lifetime: 1.5 days (top), 3 days 
(middle) and infinite (bottom). For each mode lifetime, the curves for different i and I are all normalized to have the same area. 

Finally, we can express these in terms of Dq under the assumption that the asymptotic relation (equation [T]) holds 
exactly, although in fact this is not likely to be the case: 

^i^Jven.odd = 0.04i?o (B7) 

and 

'^'^cven.odd = 0.06i?o. (B8) 



